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We study the characteristics of the tensor correlation in 4 He using a shell model type 



method. We treat the tensor force explicitly by performing a configuration-mixing calcula- 
tion in the 2p2h basis and include single-particle states up to intermediately high angular 
momenta. We adopt the Gaussian expansion method for the quantitative description of 
the spatial shrinkage of the single-particle states to optimize the tensor correlation. We are 
able to describe the full strength of the tensor correlation for 4 He in the shell model type 
method by realizing convergence. We call this model the tensor-optimized shell model. It 
is found that in 4 He, three specific 2p2h configurations are strongly coupled with the (0s) 4 
configuration due to the characteristic features of the tensor operator. 

§1. Introduction 



O 

The tensor force is an essential component of the nuclear force and plays a 
significant role in nuclear structure. For 4 He, accurate calculations using a realistic 
nucleon-nucleon interaction demonstrated that the contribution of the tensor force 
is very large, reaching —68 MeV, which is of the same order as or even larger than 
that of the central force. 1 )~ 3 * ) The resulting strong tensor correlation enhances the In- 
state probability up to about 15% in the wave function. In light nuclei up to around 
mass number 10 (A < 10), the Green's function Monte Cairo (GFMC) method 4 ) 
was developed to treat realistic interactions and is capable of describing ground 
states and a few excited states. With this method, it was demonstrated that the 
contribution of the one-pion exchange potential (OPEP) provides 70%-80% of the 
two-body attraction, in which the dominant component of OPEP is the tensor force. 

We wish to describe nuclei with larger mass numbers explicitly using realistic 
nucleon-nucleon interactions. Variational methods including GFMC of few-body 
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systems use the relative coordinates of nucleons, which are suitable to work out the 
nucleon-nucleon interaction. We call this the many-body theory with T-type basis, 
because the two interacting nucleons, whose centers of mass are connected by some 
reference coordinate, interact directly through the nucleon-nucleon interaction. In 
this case, the number of relative coordinates increase as A(A — l)/2, where A is the 
number of nucleons. Hence, the T-type basis is not suited for heavier systems, and 
calculations become increasingly difficult as A increases. By contrast, the mean- 
field framework uses as the coordinates the positions of the nucleons relative to the 
center of the nucleus. In this case, the number of coordinates is A, and hence this 
framework is suited for heavy systems. We call this method the many-body theory 
with V-type basis, since the two interacting nucleons are labeled by the coordinates 
representing the distances from the center of the nucleus, and they interact only 
indirectly through the nucleon-nucleon interaction. Hence, it is difficult to treat 
realistic interactions involving short-range repulsion and the tensor force, because 
many configurations in the V-type basis are needed to describe motion under such 
a realistic interaction. In this case, we have to invent some method to treat these 
features of realistic interactions. 

The standard method for the description of many-body systems in the V-type 
basis is the Brueckner-Hartree-Fock (BHF) method, 5 - ) ~ 7 - ) in which two-body inter- 
actions with a short-range repulsion and a strong tensor force are treated with the 
Brueckner G-matrix effective interaction under the independent pair approximation. 
Since the resulting G-matrix effective interaction is a smooth interaction, the wave 
functions in the V-type basis can treat the G-matrix within the model space. The 
BHF method or the shell model, which is based on the Brueckner theory, have some 
success in the description of many-body systems. However, there seem to be some 
essential features missing, such as the spin-orbit interaction 8 )' 9 ) and we are therefore 
forced to employ some phenomenology. 

Recently, there have been two important steps proposed for the full description 
of nuclei in the V-type basis. One is the method of Neff and Feldmeier which 
treats the unitary transformation of the radial correlation due to the short-range 
repulsion and the tensor correlation due to the tensor force separately. 10 ) The unitary 
transformation method introduces the unitary operator for the short-range repulsion 
and another unitary operator for the tensor force. It seems that the short-range 
repulsion is treated properly, because it acts at a very small distance (r < 0.5 fm), 
and therefore the terms more than the three-body operators caused by the radial 
unitary transformations can be ignored. 11 ) By contrast, the tensor force is not 
treated properly within the two-body terms of the unitary transformed Hamiltonian. 
This may be related to the fact that the range of the tensor force is intermediate 
(0.5 fm < r <1.4 fm). 

The second proposed method is that of Toki et al. 12 ) for the treatment of the pion 
exchange interaction in terms of the relativistic mean-field framework. The mean- 
field method was introduced in the relativistic mean-field approximation with a finite 
pion mean field. Because the pion is a pseudo-scalar meson, the finite pion mean field 
in the spherical ansatz naturally induces parity mixing of the intrinsic single-particle 
states. This relativistic mean-field model provides an interesting mechanism for the 
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splitting of the spin-orbit partner due to the mixing of the parity partners with the 
same total spin as Si/ 2 and pi/ 2 or p 3 / 2 and d 3 / 2 . 13 ^ The results of this work led 
us to consider employing the parity mixing single-particle states for the description 
of the tensor correlation in the non-relativistic approximation. Furthermore, it is 
straightforward to perform the parity projection and the charge projection, due to 
the fact that the tensor force is isospin dependent. 14 ) The leading term in the parity 
and charge projection consists of 2p2h configurations. 12 ) 

With the results described above, the next step was to employ a shell model 
type prescription in which the closed shell wave function is treated as the OpOh state 
and 2p2h wave functions are introduced to treat the tensor force. 15 ) When this 
method is applied to 4 He, for which the closed shell state is (0s!/ 2 ) 4 and the particle 
states of the 2p2h configurations are those in the p-orbit, the major 2p2h state is 
(0pi/2) 2 (0s 1 / 2 ) _2 ) pionic state with the CP coupling in the particle-hole picture, 12 ) 
and the spatial extension of the 0p±/ 2 wave function is almost half that of the ordinary 
harmonic oscillator wave function. 15 ) In this case, however, we still had to strengthen 
the matrix elements of the tensor force by 50% in order to provide a sufficiently large 
amount of the D-state probability of 4 He. We have to understand why we ought to 
increase the tensor matrix elements by 50%. 

In this paper, we would like to extend further the shell model type method by 
increasing the number of 2p2h configurations in order to understand the origin of 
the 50% increase of the tensor matrix element. Let us consider more carefully the 
characteristics of the tensor force. The tensor force is Vp = fT(r)S\ 2 , with the 
tensor operator S± 2 = \/24ir[Y 2 (ri 2 ), [a 1 , <r 2 ] 2 ] . Both the coupled spin (S) and the 
coupled orbital angular momentum (L) of two particles must flip with AS = AL = 2. 
Hence, the lowest relative angular momentum state for two particles is the L = 2 
state, for which the centrifugal potential is given by V cen = 6/j 2 /(mjvr 2 ), and it is 
1 GeV at 0.5 fm. This means we should describe the distance as short as about 
0.5 fm and larger by introducing higher angular momentum states. We have the 
relation /t(?") = Yli fi(ri,r 2 )Pi(fi 2 ). Here, Pi(r\ 2 ) has a significant strength around 
f ~ 0, with a width of about l/l for large / (7 > 5). We can write Pi(fi 2 ) = 
(4ir/\/2l + 1)(— y[Yi(fi), Yi(r 2 )]o. Because the relative distance to be expressed in 
terms of the high-spin state with / of the V-type coordinate is r\ x (l/l) > 0.5 fm, 
with r\ ~ r 2 ~ 2 fm (the size of 4 He), we find I < 4. In addition, there is an 
additional change of the angular momentum by 2 for the case of the tensor force. 
This means that we may be able to describe the tensor matrix elements by taking 
the single-particle states with angular momentum up to / = 4 + 2 = 6 for the 4 He 
case. Hence, we may be able to describe the tensor force in terms of the T^-type 
coordinates with intermediate size angular momenta. It remains to construct a shell 
model type prescription to handle the tensor force by using the V-type basis. We 
name this model the tensor optimized shell model. 

The tensor optimized shell model (TOSM) is described as follows. 

(1) We base it on a closed shell state, like the OpOh state, and add 2p2h states. 

(2) We include intermediate high angular momentum orbits for the particle states 
within the 2p2h excitations. 
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(3) We allow radial shrinkage of the particle states by employing the Gaussian wave 
function with the lowest node for each angular momentum state. 

(4) We use the many-body Hamiltonian with a suitably chosen two-body interaction 
that explicitly includes the tensor force. 



the particle states of finite angular momenta as large as the result converges, with 
convergence expected to occur around 1 = 6. No angular momentum states have 
nodes, except the s-state, which is constructed as a state orthogonal to the Os state 
with Gaussian wave functions with various size parameters. We do not consider the 
problem for the case of a short-range repulsion in this paper and instead use one 
obtained phenomenologically. However, once we have succeeded in describing the 
tensor force, we plan to also determine the short-range part. 

In ^21 we describe in detail the tensor-optimized shell model (TOSM), which 
treats the tensor correlation explicitly. In $3] we investigate the structure of 4 He 
by considering high angular momentum states and discuss the characteristics of 
the tensor correlation. In 21 we study the tensor correlation in detail in terms of 
the radial shrinkage with the Gaussian expansion method. In £j5l we treat 4 He by 
employing the bare tensor force of AV8' with a slight modification of the central 
interaction in the TOSM. The present study is summarized in £}6j 



In this paper, we apply the tensor optimized shell model to 4 He. Particularly, 
we would like to see if the tensor correlation is satisfactorily described by employing 
Gaussian wave functions with intermediate size angular momenta. 

2.1. Higher partial waves and variational methods 

We base 2p2h wave functions in addition to the closed shell wave function as 
the OpOh state and adopt the (Os^) 4 state for 4 He. We include 2p2h excitations 
with higher partial waves up to / ma x> the maximum orbital angular momentum of 
the excited particle orbits. We successively increase the value of Z max and watch the 
convergence of the solutions. The Is orbit is included in the cases Z max > 2. We de- 
scribe the particle states while maintaining their orthogonality to the occupied states. 
We do not include higher nodal orbits, such as 2s and lp, except for the Is orbit. 
The reason is as follows; We have confirmed that with the usual shell model type 
calculation including the higher nodal orbits of harmonic oscillator wave functions 
(HOWF) with a common length parameter, it is difficult to describe the tensor cor- 
relation satisfactorily. In fact, the tensor correlation was investigated using a higher 
nodal HOWF in a previous study, 16 ) and it was shown that the convergence of the 
solutions with respect to the principal quantum number N is slow. Furthermore, 
we found that the tensor correlation is optimized with the small size parameters of 
the Gaussian wave functions for the excited particle states. These results imply that 
the ordinary HOWF is not sufficient to satisfactorily represent the tensor correlation. 
For this reason, instead of HOWF, we adopt the Gaussian expansion technique 17 )' 18 ** 
for single-particle orbits. Each Gaussian basis function has the form of a nodeless 




state as the basic state and add 2p2h states with 



2. Tensor optimized shell model for 4 He 
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HOWF, except for Is orbit, and when we superpose a sufficient number of Gaussian 
basis functions with appropriate length parameters, we can fully optimize the radial 
component of every orbit. 

The form of the Gaussian basis function for the label a, such as Os^ and Qpi/2, 
is expressed as 

Mr,b a , m ) = Niiba.m) r l e -(r/ b «.™) 2 /2 [Y^Xi^j, (2-1) 



Ni(b a>m ) 



9 t-(2i+3) 
^ Oa,m 

T(Z + 3/2) 



(2-2) 



where m is an index that distinguish Gaussian basis functions with different values 
of the length parameter b a ^ m . We construct the following ortho-normalized single- 
particle wave function ip™ with a linear combination of Gaussian bases 

^(r) = J2 d %,™ ( t><*(r,b a ,m) for n = l,~',N a . (2-3) 

m=l 

Here, N a is a number of basis functions for a. The coefficients {<i™ m } are determined 
by solving the eigenvalue problem for the norm matrix of the non-orthogonal Gaus- 
sian basis set given in Eq. (|2-1|) , We obtain the new single-particle wave functions 
{4>a} i n Eq. (|2-3p which are distinguished by the index n for different values of the 
radial component with a given label a. 

We choose the Gaussian basis functions for the particle states to be orthogonal 
to the occupied single-particle states, which are Osi/2 in the case of 4 He. For Osi/2 
states, we simply employ one Gaussian basis function, namely, HOWF with length 
bos 1/2 ,m=i = bos- For lsi/2 states, we introduce a basis orthogonal to the 0si/ 2 states 
and possessing a length parameter b\ S)m that differs from bo s . This extended lsi/2 
state is expressed as 

<t>ls{r, bos, bu,w) = N(bis t m) {fo(b ) + f2(bos,hs,m) r } 

x e -W£>is,™) 2 /2 Yoo (r), (2-4) 

where we omit the spin wave function. Condition of coefficients /o and /2 are de- 
termined by the normalization of <f)\ s and the orthogonality to the Os state. The 
ortho-normalized Is wave functions {ipi s } are obtained from <pi s with various range 
parameters, as described above. 

Here we set the variational wave function for 4 He. Each configuration of 4 He is 
expressed as anti-symmetrized products of the four single-particle states, which are 
distinguished by a and the basis index re, as 

*p = AflrA, (2-5) 

fc=l J 

where p indicates the set re&; k = 1, • • ■ ,4}. The total wave function of 4 He with 
(J n ,T) = (0 + ,0) for spin and isospin is expressed a superposition of wave function 
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<P P as 



!P-( 4 He) = a P $p, (2-6) 



where a p represents variational coefficients. 

The variation of the energy expectation value with respect to the total wave 
function !^( 4 He) is given by 

*S = o ■ M 

which leads to the following equations: 

d{*\H-E\#) n a<^|g-^) n , 9R , 

HI ~~ U ' a ~~ U • \ Z '° 

oo a ,m Oa P 

Here, E is a Lagrange multiplier corresponding to the total energy. The parameters 
{b a , m} for the Gaussian bases appear in non-linear forms in the energy expectation 
value. We solve these two kinds of variational equations in the following steps. 
First, fixing all the length parameters b a>m , we solve the linear equation for {a p } as 
an eigenvalue problem for H with partial waves up to Z max - We thereby obtain the 
eigenvalue E, which is a function of {6 a ,m}- Next, we try various sets of the length 
parameters {b a ^ m } in order to find the solution which minimizes the total energy. In 
this wave function, we can describe the spatial shrinkage with an appropriate radial 
form, which is important for the tensor correlation. 

The analysis carried out in this paper consists of two steps. The first is the 
main analysis in which we prepare a single Gaussian basis function for each label 
oi [N a = 1 in Eq. (|2-3|) ]. In this case, we study the convergence of the energy by 
including higher partial waves and investigate the characteristics of the configuration 
mixing. As the second part of the analysis, we adopt the Gaussian expansion method 
for the excited nucleon states and examine the quantitative description of the spatial 
shrinkage of the wave function. 

Our wave function may contain an excitation of the spurious center of mass (cm) 
motion. We estimate the amount of this excitation using a following operator Hq; 

H G = T G + ±M G u 2 R G , T G = (2-9) 

l A A h 

M G = A-m N , RG = ^Y j r l , P G = V Pi , u = - r , (2-10) 

i=l i=l m N0 Qs 



where A is a mass number. When the cm. motion is constrained to the 0s state 
in the total wave function, (H G ) must be 1.5 tvjj. We estimate the amount of the 
spurious cm motion by considering the deviation of (H G ) from 1.5huj. 



Tensor Correlation in 4 He with the Tensor- Optimized Shell Model 



7 



2.2. Hamiltonian 

The Hamiltonian for 4 He is given by 

4 4 

H = J2u-T G + J2v ij , (2-11) 

i=l i<j 

with 

Vij=v C +v T +v LS + v Clmb (2 . 12) 

The effective NN interaction Vij consists of central {vfA, tensor (vfj), LS (v^ s ) and 
Coulomb (vj< lmb ) terms. We do not explicitly treat the short-range correlation and 
instead use some phenomenological central interactions. For the tensor force, we 
use one of the effective tensor forces in which the short-range part is renormalized, 
as this has been the standard choice for our study of the tensor correlation. We 
demonstrate toward the end of this paper, however, that the tensor correlation does 
not depend on the effective tensor interaction qualitatively, because of the strong 
centrifugal potential. 15 ) 

We next explain effective interactions that include the tensor force explicitly 
and are suitable for the present analysis. We use a sufficiently strong tensor force 
comparable to realistic interactions, because we would like to determine the charac- 
teristics of the tensor correlation in the structure of 4 He. In Ref.15), we examined 
two kinds of effective interaction, AK 19 )' 20 ) and GPT. 21 ) The former is constructed 
from the AV8' realistic interaction using the G-matrix theory, where the cutoff mo- 
mentum of the Q-space for the G-matrices is chosen as ItQ = 2.8 fm _1 , twice the 
usual Fermi momentum (kp = 1.4 fm _1 ). In this prescription, short-range correla- 
tions including very high momentum components (k > kg) of the tensor correlations 
are renormalized into the central term of the G-matrices. The tensor force giving 
momentum components with k < ItQ in the G-matrices survives in intermediate and 
long ranges. 19 )' 20 ) In fact, the matrix element of the tensor force given by AK is very 
close to that obtained using the bare case (AV8'), which is discussed in ^5) We also 
consider the so-called GPT interaction, 21 ) which has a mild short-range repulsion 
and was constructed to reproduce the two-nucleon properties. In GPT, renormal- 
ization of the tensor force is regarded as being based on the conventional approach, 
which gives a weak tensor force. 

We found that AK gives a large tensor contribution in 4 He, 15 ) and it seems to 
be suitable for the present model to investigate the tensor correlation. However, as 
shown below, this force gives a small radius and results in an overbinding problem for 
4 He. Contrastingly, GPT gives a good radius and a good binding energy, but a small 
tensor contribution. 15 ) It is necessary to construct other appropriate interactions 
that take into account the short-range and tensor correlations consistently for the 
present model. This is, however, beyond the scope of the present paper. In this 
study, we phenomenologically improve AK, with the central part replaced by that 
of GPT; we call this interaction GA. With this interaction, the total energy and the 
radius are improved in 4 He. 15 ) Thus, GA can be regarded as a phenomenological 
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effective interaction, and we use this interaction to study the tensor correlation for 
4 He. 

§3. Results of tensor correlation in the TOSM 

We present here the results of the configuration mixing in TOSM and the con- 
tributions of higher angular momentum states for the tensor correlation in 4 He. We 
first present the calculated results with partial waves for £ max =l to 6, with a single 
Gaussian basis (N a = 1) for each single-particle state tp a , where we optimize the 
length parameters of the Gaussian basis {b a ^ m= i} = {b a } for every state. We also 
discuss the characteristic features of the tensor correlation. 

3.1. The calculated results with intermediate high angular momentum states 

Here we present the results for the energy minimum point in each case of / max . 
The results for the set {b a } are listed in Table U and for the properties of the He 
wave functions in Table [Hi From Table HI the energy minima are found to satisfy 
ba^Os ~ 0-6 bos, which indicates spatial shrinkage of the wave function. As discussed 
below, the expectation value of the tensor force (Vr) gives the largest contribution 
around this energy minimum. 

From Table |H] and Fig. [IJ with Z max we see the convergence of the energy, the 
central force contribution (Vc), {Vr), and the D-state probability P{D). Here, we 
define P{D) as the component for which the values of the total orbital angular 
momentum and total spin are both 2 in the total wave function l^( 4 He) in Eq. (|2-6p . 
It is also found that the convergence of {Vr) is slower than that of {Vc)- This implies 
that the higher partial waves are necessary for the tensor correlation. For (Vr) and 
P(D), the I max — 4 case is almost sufficient in the present model. The excitation of 
the spurious cm. motion is small, as judged from the values of {Hq), and decreases 
with Z max as shown in Table HT1 

In Table IIIH we list values for the six configurations in the order of their proba- 
bilities for each case of Z max . For 2p2h states, excited two-particle states are shown. 

Table I. The optimized length parameters of single-particle states for 4 He in unit of fm. 



^max 


0Sl/2 


0pi/2 


0P3/2 


lSl/2 


0d 3 /2 


0d 5/2 


O/5/2 


O/7/2 


l 


1.26 


0.75 


0.69 












2 


1.19 


0.78 


0.75 


0.76 


0.69 


0.62 






3 


1.16 


0.74 


0.66 


0.73 


0.67 


0.62 


0.77 


0.66 


4 


1.16 


0.75 


0.67 


0.73 


0.67 


0.61 


0.77 


0.67 


5 


1.16 


0.76 


0.67 


0.73 


0.67 


0.61 


0.77 


0.64 


6 


1.16 


0.76 


0.67 


0.73 


0.67 


0.61 


0.77 


0.64 



^max 


O37/2 


O39/2 


Oft.9/2 


Ohu/2 


Oin/2 


0il3/2 


4 


0.70 


0.67 










5 


0.70 


0.68 


0.70 


0.67 






6 


0.70 


0.67 


0.70 


0.69 


0.71 


0.65 
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The subscripts 00, 01, 10 and 11 represent J and T, the spin and isospin for the 
two-nucleon pair, respectively. It is found that the three kinds of 2p2h configura- 
tions with (0p 1/2 )? , [(lsi/ 2 )(Od 3 / 2 )]io and [(Op 3 / 2 )(0/ 5 / 2 )]io for the particle part are 
significantly mixed. They all have the values (J, T)=(1,0), which are the same as 
those for the deuteron, and thus this two-nucleon coupling can be understood as a 
deuteron-like correlation. 15 ) 

In order to see the relation between the above three configurations and the tensor 
correlation, we expand (Vt) into two terms of matrix elements between OpOh and 
2p2h components and between 2p2h and 2p2h components as follows: 

(V T ) = ^a P a p > ($ p \V T \<P p/ } = {0p0h-2p2h} + {2p2h-2p2h}, (3-1) 

p,p' 



Table II. The properties of 4 He for each Z m ax- We list the contributions from each term in the 
Hamiltonian (in units of MeV) and Hg (in 1.5 hu>), the matter radii (i? m , in fm) and P(D) (in 
%). 



^max 


1 


2 


3 


4 


5 


6 


E 


-22.66 


-33.82 


-40.85 


-43.65 


-44.55 


-44.85 


(T> 


66.77 


82.25 


90.53 


91.83 


92.08 


92.16 


(Vc) 


-68.59 


-75.57 


-79.16 


-79.46 


-79.43 


-79.40 


(V T ) 


-21.65 


-41.17 


-52.80 


-56.52 


-57.67 


-58.06 


(Vls) 


-0.67 


-0.20 


-0.30 


-0.36 


-0.39 


-0.40 


(Vdmb) 


0.87 


0.87 


0.88 


0.87 


0.86 


0.86 


Rm 


1.34 


1.27 


1.24 


1.24 


1.24 


1.24 


P(D) 


4.66 


6.78 


7.73 


7.82 


7.85 


7.83 


<tf G }/(1.5M 


1.067 


1.048 


1.024 


1.021 


1.017 


1.016 
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Table III. The configurations of 4 He in order of decreasing probabilities for each case. 



^max 1 


/ —9 

"max 


/ —3 


OpOh 
(Opi/ 2 )?o 

(0p 3 / 2 )?0 

(Opi/ 2 )oi 

(Op 3 / 2 )oi 

(Opi/ 2 )(Op3/ 2 )io 


91.61 
6.16 
1.25 
0.55 
0.37 
0.07 


OpOh 

(Opi/ 2 )io 
(ls 1/2 )(0d 3/2 )io 
(0d 3/2 )? 
(0p 3/2 )? 
(0d 5/2 )? 


89.38 
4.31 
3.00 
0.79 
0.75 
0.41 


OpOh 

(0pi/ 2 )? 
(0p 3 / 2 )(0/ 5/2 )i 
(ls 1/2 )(Od 3/2 )io 
(0p 3 / 2 )?o 
(0d 3/2 )? 


88.53 
3.82 
2.31 
2.03 
0.57 
0.54 




^max ' — 4 


^max — 5 


^max — 6 


OpOh 

(Opi/ 2 )io 
(ls 1/2 )(Od 3/2 )io 
(OP3/2) (0/5/2)10 
(Od 5/2 )(Off 7 / 2 )io 

(Od 3 / 2 )?o 


88.47 
3.62 
2.03 
1.92 
0.61 
0.54 


OpOh 

(0pi/ 2 )u> 
(ls 1/2 )(0d 3/2 )i 
(0p 3/2 )(0/ 5/2 )i O 

(0d 3 / 2 )io 
(0d 5/2 )(0g 7/2 )i 


88.45 
3.58 
2.03 
1.90 
0.54 
0.51 


OpOh 
(Opi/ 2 )io 
(lsi/ 2 )(Od 3/2 )io 
(OP3/2) (0/5/2)10 

(0d 3/2 )? 
(Od 5 / 2 )(Off7/ 2 )io 


88.48 
3.56 
2.03 
1.90 
0.54 
0.52 



Table IV. Expansion of (Vr) into the two terms in Eqs. (|3T[) in units of MeV. 



^max 


1 


2 


3 


4 


5 


6 


0p0h-2p2h 


-21.02 


-42.03 


-54.74 


-59.25 


-60.70 


-61.22 


2p2h-2p2h 


-0.63 


0.86 


1.94 


2.73 


3.03 


3.16 



{0p0h-2 P 2h} = a o P oh a p ((^opohlVrl*,) + (*pI^t|*qpoa)) , (3-2) 

pg2p2/i 

{2p2h-2p2h} = ^2 a p a p , {<P p \V T \$ p i) ■ (3-3) 

P , P '£2p2h 

We list the contributions of these two terms in Table IIVL It is found that the first 
term {0p0h-2p2h} is dominant, and this result is consistent with that obtained in 
previous studies. 1 ) We furthermore list each component in {0p0h-2p2h} for the 
Im&x = 6 case in Table El Three configurations with large probabilities in Table IIIII 
give large contributions to (Vr), about 70% of the total value as shown in Table IVl 
This means that three configurations are essential for the tensor correlation in 4 He. 
When we increase Z max from 3 to 6, these configurations remain strongly mixed. This 
result indicates that the £ m ax = 3 case is sufficient to describe the characteristics of 
the configuration mixing for the tensor correlation in 4 He. 

3.2. Special feature and the strong shrinkage of the tensor correlation 

We next discuss the reason that three configurations are related to the tensor 
correlation. We depict the coupling scheme generated by the tensor force for the 
OpOh and these 2p2h configurations in Fig. [21 where and {sj} (i = 1, ■ • ■ ,4) are 
the orbital angular momenta and the directions of the intrinsic spins for each nucleon 
in 4 He, respectively. The OpOh configuration of (0s) 4 is expressed as 1/V2 {(Osi/2)oi 
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Table V. Contributions of each 0p0h-2p2h coupling to (Vt) in MeV for the I max — O CEtSC. 



2p 


(&0p0h\V T \$2p2h) 


(Opi/ 2 )io 


-14.47 


(ls 1/2 )(Od 3/2 )io 


-15.01 


(QP3/2) (0/5/2)10 


-12.00 


(Od 3/2 )?o 


-3.67 


(0d 5 / 2 )(037/ 2 )io 


-5.17 



(Osi/2)oi + (O s i/2)io(O s i/2)io}- When the tensor force acts on (0s 1 / 2 )f in the OpOh 
configuration, it changes both the orbital angular momentum (L) and intrinsic spin 
(S) for excited particle states by 2 {AL = 2 and AS = 2). The directions of AL and 
AS must be anti-parallel, since the rank of the tensor operator 5*12 oc [a 1 , <x 2 ] 2 ]o 
is zero. The direction of S in the particle states of the 2p2h configuration is also 
opposite to that of the original (0si/ 2 ) 2 n component in the OpOh case. If the L = 2 
component in the 2p2h configuration consists of Z3 = 1 4 = 1 for orbital angular 
momenta of two excited nucleons, as shown in Fig. [2j the directions of the orbital 
angular momentum and the spin for each excited single-particle state will be op- 
posite, and hence this state becomes 0pi/ 2 and (Opi/ 2 )io ^ s formed. The situation 
discussed here can also be understood as the 0~ coupling between 0si/ 2 and 0pi/ 2 
states in the particle- hole picture. 12 )~ 14 ) If the L = 2 component consists of I3 = 
and Z4 = 2, the directions of Z4 and the spin for the single-particle state are also 
opposite, and this state becomes 0^3/2- This component makes [(lsi/ 2 )(Od3/2)]io as 
shown in Fig. [2j If I3 = 1 and U = 3, [(0p 3 / 2 ) (0/5/ 2 )]io is obtained in a similar 
way. Higher configurations, having a set of larger orbital angular momenta of I3 
and Z4 possess small amplitudes, due to their large kinetic energies. For this reason, 



*He(0 + ) 



2p2h : [(Os 1/2 )io(Opi/2)io]oo- 



| 1= | 2 =0 L=2A 

4 s=i T 

[(Osi/ 2 )fo(Os 1/2 )f ]oo S ^ S4 Hi 



OpOh : (0s 1/2 ; 

ZD _ 

[ 1= 1 2= | 3= | 4= 

Slj|l S 2 S3ftt s 4 

. 4 = 1 + J2 = 1 + . 
V 

J =J 1 ®J 2 =0 



13=14=1 



jp = r 



selectivity of tensor operator 
AL=2 AS=2 with AL®AS=0 




2p2h : 

( [(OSi/2)?o[(1Sl/2)(Od 3/ 2)]lo]oo, 

| 1= | 2 =0 L=2 , 
S=1 

si|||s 2 s 3 ijis 4 | 3= o| U=2 
J?=1 + jj[=1 + 



2p2h:[(Os 1/2 )? [(Op3 /2 )(Of 5/2 )] 1 J]oo 



h=l 2 =o s=1 L=2- 

si4|s 2 s 3 i|is 4 l 3 =l||| 
J?=1 + J^=1 + 



Fig. 2. Schematic depiction of the favored coupling through the tensor force for He. 
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2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 2 4 6 8 10 12 14 16 18 
N=2n+I [hco] N=2n+I [hco] N=2n+I [Moo] 



Fig. 3. Accumulative probabilities of the single-particle states for 4 He in terms of HOWF for the 
fmax- 4 case. 



100 




0.5 _ 1.0 1.5 

b/b 0s 

Fig. 4. The properties of 4 He in the Z ma x=4 case as a function of b/bo s - P2 P 2h(JT) are the proba- 
bilities of 2p2h components with JT being the spin-isospin pair of the excited two nucleons. 

the strong mixing of three 2p2h configurations in 4 He can be understood from the 
selectivity of the operator of the tensor force. 

Here we discuss the relation of the spatial shrinkage of the single-particle states 
and the higher shell effect. For this purpose, we expand the shrunk single-particle 
wave function in terms of HOWF in the same length parameter of the 0s state. 
We calculate accumulative probability defined as the summation of the expansion 
probabilities up to the quanta N = 2n + 1. The results are shown in Fig. [3] for the 
Imax = 4 case. A higher shell effect is seen, as 2p2h configurations require at least 
\AJ\w excitations for the shrunk states. 

In Fig. HI we plot the length parameter dependence of the solutions of 4 He in 
order to see the effect of the spatial shrinkage more clearly, where length parameters, 
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except for bo s , are changed as a common value b. We plot the solutions as a function 
of the length parameter ratio b/bo s in the Z max = 4 case, where bo s is fixed to 1.16 fm. 
The energy minimum is found to be near b = 0.6 bo s , where (Vr) gives the largest 
contribution and the kinetic energy takes its maximum value. In this case, (Vr) at 
the energy minimum is more than twice as large as in the ordinary shell model case 
with b/bo. s =l. For the 2p2h configurations, the (J,T) = (1,0) pair is strongly mixed 
around the energy minimum. From these results, the tensor force is incorporated 
with a small length parameter for the excited nucleon states. 

As shown in Table HI the obtained wave functions give matter radii that are 
smaller than the experimental value (1.48 fm) by 0.2 fm. This is due to the fact that 
in addition to the overbinding of the system, the length parameters of all single- 
particle states are smaller than in the case of the single configuration of (0s) for 
4 He with bo s =lA fm, which reproduces the observed matter radii without the tensor 
correlation. 

§4. Detailed analysis of the radial wave functions 

In this section, we improve the radial wave functions of the particle states with 
the Gaussian expansion technique in order to obtain a quantitative description of 
the tensor correlation. We adopt the Gaussian expansion method 17 )' 18 - 1 for every 
excited particle state. This technique has been shown to possess wide applicability 
for describing short-range correlations in few-body calculations and weakly bound 
states in neutron-rich nuclei. 

We now present the results in the l m ax=2 case. Here, we fix bo s to 1.4 fm to 
obtain the appropriate matter radius of 4 He in the present wave function, as the 
results for the radius listed in Table [TT] are too small. We use the same number of 
basis functions, N a , for every a of the particle states. The set of length parameters 
for the Gaussian basis functions in Eq. (|2-ip are listed in Table |VH 

From Fig. we confirm the convergence of the solutions with respect to N a . In 
particular, the physical quantities almost converge within three or four bases, where 
the bases with length parameters smaller than bo s are favored. It is also found that 
the Gaussian expansion gives a larger effect on (Vr) than (Vc)- In particular, P{D) 
attains a value nearly twice as large as that in the single Gaussian case (N a = 1). 
The kinetic energy also increases about 20% in the Gaussian expansion in comparison 
with the N a = 1 case. The matter radius converges to 1.39 fm. We thus find that the 
Gaussian expansion succeeds in describing the radial wave functions for the excited 
particle states and increases the tensor correlation for 4 He. In order to estimate this 
effect, we calculate 4 He with a single Gaussian basis function. Doing so, we find 
that if we enhance the matrix elements of the tensor force by 50%, the result of the 
Gaussian expansion is simulated. This result implies that the reason for using the 
above enhanced tensor matrix elements in previous study 15 ) is confirmed from the 
improvement of the spatial description of the particle states. 
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Table VI. The set of length parameters for the 
Gaussian basis functions. 



N a 


Set Of b a ,m 


(to 


= 1, 




,N a ) 


1 


0.7 












2 


0.6 0.8 












3 


0.6 0.8 


1.0 










4 


0.6 0.8 


1.0 


1.4 








5 


0.6 0.8 


1.0 


1.4 


2.0 






6 


0.4 0.6 


0.8 


1.0 


1.4 


2.0 




7 


0.4 0.6 


0.8 


1.0 


1.4 


2.0 


2.5 



> 

CD 



CD 

X 



O) 
CD 

c 
LU 




Q 
Q. 



1 2 3 4 5 6 7 

Number of Gaussian basis functions 

Fig. 5. The properties of 4 He obtained with 
the Gaussian expansion method. The circle 
indicates P{D) with a single Gaussian ba- 
sis function realized by increasing the ma- 
trix elements of the tensor force by 50%. 



§5. Properties of 4 He with the tensor component of the 

nucleon-nucleon force 

In this section we present the results obtained with TOSM for 4 He using the 
tensor force of the nucleon-nucleon interaction together with the central interaction 
modified by accounting for the short-range repulsion. For this purpose, we first 
compare the matrix elements of the tensor forces obtained in various methods. 

5.1. The matrix elements of various tensor forces 

We compare the radial matrix elements of the tensor force Vt = /t(0 ■ <St2 
for several interactions: the bare nucleon-nucleon interaction, AV8', 22 ) the Bonn 
potential with it + p, 23 ^ and the G-matrix tensor interaction, AK, which is based on 
AV8'. In Fig. [7J we show the radial dependence /r(^) of each tensor force for the 
triplet-even channel. We see that the short-range behavior depends very strongly on 
the choice of the nucleon-nucleon interaction. Contrastingly, the intermediate- and 
long-range behavior depends only weakly on the choice of the potential. 

We now compare the radial integrals for three kinds of tensor forces for the S-D 
coupling in the relative motion defined by the following equations: 

/>oo 

/ SD = / dr MsdW, M SD (r) = r 2 cp 0s (r, b 0s ) ■ f T (r) ■ cP od (r, b od ). (5-1) 
Jo 

This matrix element is essential to explicitly describe the tensor correlation. Here, 
4>o s and (pod are the radial components of the Gaussian wave functions for the s and d 
waves with the length parameters 60s and &od> respectively. We choose bo s = 1.4- \[2 
fm for the 0s wave function, where the factor \f2 comes from the relative coordinate 
transformation of HOWF. We consider two length parameters, &od = and 
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bod = bo s , for the Od wave function. The smaller length parameter bod, corresponds to 
the case of the largest matrix element of the tensor force, as discussed in the previous 
sections and the previous study, 15 ) and the larger length parameter corresponds to 
the case of the tensor force in the standard shell model. 

In the left panel of Fig. [7J we depict the shorter range case for the radial inte- 
grands Msd(^) in the three cases, which are seen to be very similar. The intermediate 
and long-range behavior is important for the matrix elements of the tensor force. The 
large difference in the short-range part is washed out due to the AL = 2 transition 
of the tensor force for the relative orbital angular momentum. Due to the centrifugal 
barrier of the D state, the contribution from short-range part becomes small in the 
tensor-force matrix elements. This means that the coupling between the short-range 
correlation and the tensor correlation is weak, and thus we can describe the tensor 
correlation explicitly in the shell model type method while we can renormalize the 
short-range correlation into the central force independently from the tensor force. In 
the right panel of Fig. [7J we also depict the longer range case. Here, again the radial 
integrands are very similar. Further, we list the matrix elements for each tensor 
force, Isd, in Table IVTil We find that the three matrix elements are very similar, 
with differences on the order of 10%. We also see that the tensor matrix element 
for the shorter length parameter is much larger than that for for the larger length 
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Table VII. Radial matrix elements of the tensor force 7sd for the triplet-even channel in units of 
MeV. 





AV8' 


AK 


7T + p (Bonn) 


bod — bo s /2 


-14.98 


-15.04 


-13.30 


bod = bos 


-3.83 


-3.85 


-3.82 



parameter, which corresponds to the case of the standard shell model. 

5.2. The properties of 4 He in TOSM 

We have studied the full treatment of the tensor force in the shell model type 
method by employing a sufficiently large number of angular momenta and radial wave 
functions using the T^-type coordinates. We have used the effective nucleon-nucleon 
interaction, GA. We found good convergence in the case of intermediate angular 
momenta with considering the several values of the range parameters for the tensor 
force. Because this is the first time that the tensor force has been treated fully, while 
the short-range repulsion has to be treated somehow at present, the binding energy 
of 4 He that we obtained is too large, and also, the size of 4 He is too small. Here, we 
would like to adjust the parameters of the central interaction so that realistic values 
for the properties of 4 He are obtained. In the previous subsection, in addition, we 
showed that we are able to treat the bare tensor force obtained from nucleon-nucleon 
scattering and the matrix elements of tensor forces are found to be very similar, due 
to the large centrifugal potential. Considering these points, we employ the bare 
tensor force for the study of 4 He and use a phenomenological central interaction in 
TOSM. 

We adopt the tensor force and the LS force of AV8'. For the central interaction, 
we begin with the central force of the GPT interaction, which is obtained by fitting 
the observed two-nucleon properties. We then adjust the central part of GPT to 
improve the large binding energy, as explained in Table HH This is carried out by 
reducing the strength of the second range (attractive part) of the GPT's central force 
by 13% in order to fit the observed binding energy (28.3 MeV) for 4 He. 

For the 4 He wave function, we choose / m ax=5 for high angular momentum states 
and use four Gaussian basis functions for every particle state, whose length param- 
eters are 0.6, 0.8, 1.0 and 1.4 fm, respectively. These choices are sufficient in the 
present calculation. We variationally determine the length parameter of the 0s state 
to be 1.42 fm. 

The results for 4 He are listed in Table IVIIII For the mixing probabilities, we 
sum the probabilities belonging to the same configurations having the same spin and 
isospin pair with different radial components obtained in the Gaussian expansion. 
It is found that (Vt) is approximately —60 MeV and P{D) reaches approximately 
9%. These results imply that the obtained wave function satisfactorily represents the 
characteristics of the tensor correlations. The excitation of the spurious cm motion 
is weak, and its size is estimated to be approximately 0.4 MeV. The three 2p2h 
configurations are strongly mixed, and the (Op 1 / 2 )io component is most favored. 
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Table VIII. The properties of 4 He in the Z max =5 case with the use of the bare tensor force. We list 
the contributions from each term in the Hamiltonian (in units of MeV) and the matter radius 
(R m , in fm) and the probabilities (in %) of each configuration labeling with the excited particle 
states. 



E 


-27.89 


(T) 


75.85 


(Vc) 


-44.57 


(V T ) 


-58.97 


(Vls) 


-0.98 


(Vcimb) 


0.78 


(tf G )/(1.5M 


1.011 


Rm 


1.47 



P(D) 


9.13 


OpOh 


85.02 


(Opi/ 2 )io 


3.40 


(ls 1/2 )(Od 3/2 )io 


2.35 


(PP3/2) (O/5/2) 10 


2.21 


(0pi/ 2 )(0p3/ 2 )l0 


1.61 


(0rf 5 / 2 ) (057/2)10 


0.83 


(Od 3/2 )?o 


0.60 


(0p 3 / 2 )l0 


0.53 



§6. Summary 

We have studied the tensor correlation of 4 He using a shell model type method 
with a V-type coordinate. We have developed the tensor optimized shell model 
(TOSM) to fully treat the tensor correlation. We have applied the TOSM to 4 He, 
where we have employed the (0s) 4 basis and included up to the two-particle two-hole 
(2p2h) excitations with particle states up to intermediate angular momenta. Particle 
states are expressed in terms of Gaussian wave functions for each partial wave, while 
the length parameters are allowed to vary so as to minimize the total energy. For the 
study of the tensor correlation, we first considered the G-matrix interaction based 
on the AV8' interaction 19 )' 20 ) and replace the central part of the interaction by those 
of the GPT. 21 ) 

We first calculated 4 He by including the intermediate high angular momentum 
states up to Z max successively, while employing only one Gaussian basis function for 
each orbit. We found that the spatial size of the particle states are significantly 
shrunk in order to optimize the tensor correlation for He. The Gaussian length 
parameters of high partial waves are found to be about 60 % of that of the 0s 
state. The energy gain from the tensor force greatly as a function of l max and 
exhibits convergence near Z max = 4. We found that three particular 2p2h states 
admix to the (0s) 4 component due to the characteristic features of the tensor force. 
We then expressed the radial wave function in terms of the superposition of Gaussian 
basis functions with several length parameters. We found that approximately three 
Gaussian basis functions are sufficient to express the radial wave functions, and the 
energy gain resulting from the use of various Gaussian wave functions is about 50% in 
comparison with the case of one Gaussian wave function. Hence, we have succeeded 
in developing a shell model type model based on the V-type coordinate, TOSM, 
which is capable of treating a strong tensor force. 

We then studied 4 He with TOSM using the bare tensor force AV8' and the short- 
range modified central interaction GPT. We were able to determine the ground state 
properties of 4 He with a modification of the intermediate range of GPT central in- 
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teraction and thereby reproduced the binding energy. The expectation value of the 
tensor force is very large, about —60 MeV, and the D-state probability is approxi- 
mately 10%. 

We developed the TOSM, with which we are able to treat the bare tensor force, 
with the hope of investigating the tensor correlations in heavy nuclei. As the next 
step we should develop a method to incorporate the short-range repulsion so that 
we are able to treat intermediate and heavy nuclei using the nucleon-nucleon inter- 
action. We plan also to study the effects of the tensor correlation on the structures 
of the nuclei in the neighborhood of 4 He. Particularly interesting are the scattering 
phenomena of 4 He+neutron system (A = 5) and the structures of A = 6 systems, 
6 He, 6 Li and 6 Be. 
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